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Barrierless procedure for substitutionally doping graphene sheets with boron atoms: 

ab initio calculations 
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Using ab initio methods, we propose a simple and effective way to substitutionally dope graphene 
sheets with Boron. The method consists of selectively exposing each side of the graphene sheet 
to different elements. We first expose one side of the membrane to Boron, while the other side is 
exposed to Nitrogen. Proceeding this way, the B atoms will be spontaneously incorporated into the 
graphene membrane, without any activation barrier. In a second step, the system should be exposed 
to a H-rich environment, that will remove the CN radical from the layer and form HCN, leading to 
a perfect substitutional doping. 
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Graphene recently became one of the most studied ma- 
terials in the literature due to its potential applications 
for novel devices, hydrogen storage and due to its unique 
physical and chemical propertiesii^i^i^i^. Nowadays, sim- 
ple graphene-based devices are already feasible, showing 
that the initial concepts proposed for graphene can be 
manufactured^-. On the other hand, in order for func- 
tional devices to be developed, a more precise ability to 
manipulate graphene has to be reached. In order to do 
that, new functionalization processes have to be devel- 
oped. A very effective way to achieve this functionaliza- 
tion level is through doping. As most materials have to 
be doped to build functional devices with them, it is rea- 
sonable to argue that graphene also should be doped in 
order to build functional devices with it. 

Doping of nanostructures is known to be a very difficult 
taskli^. This difficulty conies from a combination of ener- 
getic and kinetic effects. The same happens for graphene, 
with few successful cases of substitutional doping. It is 
known that defects can be created by ion-irradiation^ and 
impurities can be incorporated by pyrolitic^*' processes. 
High-Temperture pyrolitic processes, for example, are ad- 
equate for applications such as hydrogen storage, but 
are too dirty for electronic purposes. A better-controlled 
method is needed for that. Doping can also be induced 
by the adsorption of molecules on top of graphenei^, but 
precise atom substitution by foreign impurities is still a 
challenge. Here we propose a procedure to substitution- 
ally dope graphene sheets, ribbons and possibly carbon 
nanotubes. The process described here does not have any 
activation barrier, i.e., it occurs spontaneously once the 
reagents are placed in different sides of graphene. This 
system is a potential candidate for high-performance Hy- 
drogen storage materia l^°i^^ , nanosensors and spin filter 
devices^- . 

Our methodology has two steps: (i) in the first step, 
one side of the graphene membrane should be exposed 
to a B-rich environment, while in the other side a N-rich 
environment should be created. This combination will 
make the system energetically unstable, and the conse- 



quence is that a C atom will be pushed out from the 
sheet, and on its place a B atom will be incorporated; 
(m) the second step is the removal of the C-N radical 
that is bonded to the B atom. In order to do that, the 
system should be inserted into an H-rich medium, or an 
acid medium. To simulate this, we insert a H atom in the 
N-rich side of the membrane. The H atom will bind to 
the C-N, removing it from the surface and forming HCN, 
completing the B-doping process. Once again, this is a 
barrierless process. 

Although complex from an experimental point of view, 
according to recent results, our approach should be fea- 
sible. Recent measurementa^^ show that macroscopic 
graphene sheets can be created and precisely placed 
on top of a support grid for transimission electron mi- 
croscopy. This is a possible experimental apparatus nec- 
essary for the creation of the membrane in order for our 
doping process to be realized. This support could be 
used as a separator for the two environments. Other re- 
cent experiments show that it is possible to apply a pres- 
sure difference across a graphene membrane, and that 
it should be impermeable to standard gases, providing 
a unique separation barrier between two distinct regions 
that is only one atom thick^^ . Our model also shows that, 
besides separating two distinct regions, depending on the 
configuration of each medium, atom substitution should 
occur. 

To study this process, we performed spin-polarized 
total energy calculations based on the density func- 
tional theor3*i^ within the generalized gradient 
approximation^^. We also used norm conserving 
pseudopotentials^^, as implemented in the SIESTA 
code^^. We employ a double zeta basis function with 
polarization orbitals^S^ a confining energy shift of 
0.06 eV^ and a mesh cutoff energy of 250 Ry for the 
grid integration. The Brillouin zone was sampled by 
using a Monkorst-Pack^^ scheme with a (15 x 1 x 15) 
k-point sampling. Smaller k-point samplings leads 
to an incorrect description of the density of states of 
graphene, leading to possibly incorrect results. Periodic 
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boundary conditions were used in the simulations, and 
the graphene sheet was described by a (3 x 4) supercell 
with 48 carbon atoms. A vacuum region of 15 A was 
used in order to separate the layer and its images in 
the direction perpendicular to the plane. All atoms of 
the system were fully relaxed within a force convergence 
criterion of 25 me v/A. 
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FIG. 1: (Color online) Schematic representations of the opti- 
mized geometries for the adsorption of B (a, 6, c ) and N (d, 
e, f) on graphene. Figures a and d are associated with the 
hep site; b and e represent the top site while c and / stands 
for the brigde site. The clearer atoms in e and / indicate the 
carbon atoms pulled out the graphene plane. 

We first analyze the stability of a single B and a single 
N atom at the surface of a graphene layer. We studied 
three different local minima for each dopant: on top of a 
C atom; on a bridge site, in between two C atoms; and 
at the center of the hexagon (hep). For both dopants, 
B and N, the most stable adsorption position is at the 
bridge site. Fig. [T] shows the possible locations and the 
relaxed configurations for B and N. From a structural 
point of view, the adsorption of a B atom does not cause 
any significant distortion on the graphene sheet, while 
the N atom pulls two of the C atoms out of the plane. 
The environment around these C atoms becomes sp'^- 
like. This can be clearly observed in Fig. [T^ and [Tf. 
The optimized structural parameters for these geometries 
and the energy difference for each site and each dopant, 
relative to most stable conformation is shown in Table [H 

We proceed to the analysis when the graphene mem- 
brane is exposed to two different environments with one 
type of atom in each side. We expect that, when the 
graphene membrane is exposed to different mediums in 
each side, the atoms will tend to be localized at the most 
stable position for each atom, as described above. As 
for both atoms the most stable site is the bridge site, we 
placed one B atom in one side and one N atom on the 
other side of the membrane bridging theses carbon atoms. 
We then let the system relax to its equilibrium structure, 
where all the forces are minimized. The N atom binds 



TABLE I: Variation of the energy (AE) for different adsorp- 
tion sites and the distance between the dopant (D) and the 
closest carbon atom (do-c) for each site. All energies are in 
eV and all distances are in A. 



Dopant 


site AE 


do-c 




top 0.06 


1.76 


B 


hep 0.13 


2.24 




bridge 0.00 


1.86 




top 0.73 


1.50 


N 


hep 0.69 


3.30 




bridge 0.00 


1.49 



to the layer, distorting it as shown before. The bonds 
between the C atom that is pulled out of the surface 
and its neighboring C atoms will be elongated and con- 
sequently weakened. The B atom represents an exception 
to the octet rule^^, and it will be energetically more sta- 
ble when bonded to three atoms. As the C-C bonds in 
the graphene membrane are weakened, the B atom will 
be pulled towards the surface, binding strongly to these 
three C atoms and replacing the original C atom from the 
membrane. Besides the energetic gain due to the three 
bonds that the B atom makes, another positive factor is 
related to size effects: B has an atomic radius compara- 
ble to C, such that there will not have much stress in the 
lattice due to the B incorporation. This mechanism has 
no activation barrier and is schematically shown in Fig. 

m 




FIG. 2: (Color online) Representative geometries, along the 
structural optimization, for the atom substitution process (a- 
d) and for the radical removal (e, _/). Larger, green spheres 
represent B, while smaller blue spheres N. C and H are grey 
and white respectively. 

In the process described before, we were able to sub- 
stitute a C atom by a B atom in the graphene layer, 
but there is still a C-N radical bonded to the dopant 
(Fig. [2li). The physical properties of this system will 
be different than a standard doping; the radical should 
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be removed in order for the doping process to be ac- 
complished. To achieve this, we insert H atoms into the 
system. In our calculation, we have inserted an H atom 
nearby the C-N radical, as shown in Fig. The H atom 
will strongly bind to the C atom, forming HCN, which 
is spontaneously removed. The final result is a perfect 
doping of the graphene layer, as shown in Fig. [2f. As 
in the first part of the process, this mechanism does not 
have any activation barrier. We also tested the H atom 
connected with the N atom, however in this situation the 
radical was not removed. The energy difference between 
these two structures (I. H connected with the C atom of 
the radical and 2. H attached to the N atom of radical) 
is O.Q6eV being the first situation energetically favorable. 
The radical removal reaction can be described as: 

Graphene BCN + H — > Graphene b + HCN, 
where {GraphenescN) stands for Boron doped graphene 
coupled to a C-N radical plus a H atom, creating Boron- 
doped graphene [Graphenes) plus a free HCN. The re- 
action is exothermic, liberating an energy of 2.7 eV. 




FIG. 3: (Color online) Total charge densities for different 
steps of the doping process. The scale is in e/A. 

The doping mechanism can be better understood by 
looking the variation of the strength of each bond during 
the process. In Figure [3] we show a total charge density 
contour plot for different stages of the doping process. 
The red color in the scale bar indicates low densities and 
the blue color indicates high charge densities. The strong 
coupling between the N atom and graphene is clear in 
Fig. [3^. In Fig[5]3 we observe that the B atom is more 
connected with the right C coupled with N. In Fig. [3J; we 
present the plot for the last stage of the B incorporation 
and we can observe that the C-N bond is stronger than 
the B-C bond, this being part of the reason why this 
radical is removed when H is inserted. The B atom makes 



three bonds with the C atoms of the graphene lattice 
and the bond between the B and the carbon atom that 
is pushed out from the sheet is mediated by a weak lone 
pair shared by B and C atoms. 
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FIG. 4: (Color online) Total density of states for different 
doping stages. The graphs are aligned by the peak located at 
-lleV. The vertical lines indicate the Fermi energy. The inset 
is a plot of the local charge density, in the (B-C-N) structure, 
around the level created by the doping. 

In Figure |4] we present the density of states (DOS) 
of graphene (pristine). Boron doped graphene attached 
with the C-N radical (B-C-N) and Boron doped graphene 
with free HCN (B-I-HCN). For B-C-N we observe, above 
the Dirac point, the existence of an unoccupied level. 
This level, shown in the inset of Figdl is originated from 
the orbital's re-hybridization due to Boron incorporation, 
and has major contributions from the carbon atoms sur- 
rounding the B atom, as well as contributions from the 
C and N present on the radical. When the H atom is 
inserted, this level is filled by the H's electron, and the 
radical will detach from the sheet. The DOS of the fi- 
nal structure (B-I-HCN) is also presented, and we can see 
that the impurity level has disappeared. Also, the Fermi 
energy shifts down with respect to the Dirac point (where 
the DOS is zero), indicating that Boron acts like a p-type 
dopant (acceptor), increasing the DOS in the Fermi level 
for this doping concentration. 

We also analyzed in more detail the structural and 
electronic properties for the final geometry of our pro- 
cess(Fig 2f), i.e, a substitutional B atom with a HCN 
molecule nearby. This is important since, in some cases, 
molecules adsorbed to graphene can also induce dop- 
ing grapheneAi. We obtained a B-C bond length of 
1.52 A which is in good agreement with an experimental 
STM measurement of 1.59 A for the B-C bond length^'^. 
This result is also in agreement with previous theoretical 
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calculations2a, where the B-C bond length is found to be 
1.53 A. The nearest-neighbor C-C bonds are compressed 
by 0.015 A when compared to undoped graphene. This 
clearly shows that the distortions are strongly localized 
around the B site. 




FIG. 5: (Color online)DOS (upper panel) and PDOS (lower 
panels) for the final structure presented in Fig 2f. The insets 
are local charge densities associated with the energy levels. 
The Fermi energy is at zero. Because isolated levels corre- 
spond to delta fiinctions in the DOS, for display purposes, we 
have artificially broadened the DOS by 0.1 eV. 



For the final structure. Fig [5] shows the decomposi- 
tion of the total DOS. It is composed of discrete HCN 
levels (lower panel), and the doped graphene levels (mid- 
dle panel). The local charge densities for these levels is 
presented on the insets, and we can clearly observe that 
HCN just weakly binds to graphene, and that the doped 
structure represents a real p-type doping of graphene. 



The mechanism described here is not solely valid for 
the exact configuration described above. It also happens 
when the B and the N atoms are located at top sites. 
More importantly, test calculations also showed that the 
B atom can be incorporated when the N atom is substi- 
tuted by P, O, CH and NH, showing that our proposition 
is quite general and robust. 



In conclusion, performing ab initio electronic structure 
calculations we have proposed an effective way to substi- 
tutionally dope graphene. This can be done by using 
graphene as a membrane to separate two different medi- 
ums. By this process, in principle, small atoms such as B 
could be readily incorporated through a barrierless pro- 
cess when the other side of the membrane is exposed to 
an environment that structurally distorts graphene. The 
methodology described here should also be valid for dop- 
ing other Carbon nanostructures such as wide nanorib- 
bons and carbon nanotubes. 
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